Inverse photoemission in strongly correlated electron systems 
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Based on exact results for small clusters of t— J model we point out the existence of several 
distinct channels in the inverse photoemission (IPES) spectrum. Hole-like quasiparticles can either 
be annihilated completely, or leave behind a variable number of spin excitations, which formed 
the dressing cloud of the annihilated hole. In the physical parameter regime the latter processes 
carry the bulk of IPES weight and although the Fermi surface takes the form of hole pockets, the 
distribution of spectal weight including the 'spin excitation bands' is reminiscent of free electrons. 
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A frequently used concept in many-body physics is 
the 'dressing' of a particle by an excitation cloud. Sud- 
den annihilation of the particle then may leave behind 
its dressing cloud, and these 'released' excitations may 
carry excess momentum and energy. Such 'shake-off' pro- 
cesses represent a channel for particle annihilation of pure 
many-body origin, and one may expect that their relative 
importance depends on the 'degree of dressing', as mea- 
sured by the quasiparticle weight Z. To see this, let us 
consider the change of the single particle spectral func- 
tion upon removing an (infinitesimal) fraction e of elec- 
trons/spin direction. In a noninteracting system thereby 
the topmost Ne poles of the photoemission (PES) spec- 
trum cross to the inverse photoemission (IPES) spectrum 
(N denotes the system size). In a Fermi liquid the same 
will happen; depending on the magnitude of Z, however, 
the corresponding decrease of the integrated PES weight 
of Ne ■ Z may be substantially less than the change of iVe 
required to maintain the correct electron count (the inte- 
grated PES weight equals the number of electrons) . Spec- 
tral weight corresponding to Ne(l—Z) electrons therefore 
must shift from PES to IPES at momenta and energies 
off the Fermi surface, and a natural mechanism to ac- 
complish this would be the 'shake-off processes' adressed 
above. Values of Zm0.2 — 0.5, as found in the t—J and 
Hubbard model near half-filling Q, thus suggest that in 
this case in fact the bulk of IPES weight is carried by 
shake-off processes. In this manuscript we present diag- 
onalization results which substantiate this general argu- 
ment. The emerging picture of the single particle spectral 
function and Fermi surface not only resolves the appar- 
ent contradiction between mounting evidence |^|| for a 
hole pocket-like Fermi surface in these models on one 
hand and the roughly free electron-like weight distribu- 
tion in their spectral function Q on the other; it is also 
reasonably consistent with experiments on cuprate su- 
perconductors. 

For the standard t—J model [|J we study the electron 
addition spectrum 

A+(k,u>) = ]T \{p n \cl A \0 n+1 )\ 2 S(cj -<) 



where \v n ) denotes the v eigenstate with n holes (in 
particular u=0 implies the ground state (GS)) and w" 
its excitation energy relative to the n-hole GS. The oper- 
ator ct is the Fourier transform of cj a —c\ Tii _ a ). 

where cj denotes the Fermion creation operator. For 
small clusters A + (k,oj) can be evaluated exactly by the 
Lanczos algorithm, the results presented below have been 
obtained in the standard 18-site cluster 0, results for the 
16 site cluster are consistent with these. 
As a first step we consider the dependence of A+(k,oj) 
on J/t and hole number rih- A trivial dependence 
on n/i is due the requirement that the k and ui inte- 
grated IPES weight equals rih- Figure 1 therefore shows 
the IPES spectra divided by rih for the 1 and 2 hole 
ground states and various J/t. The single hole GS 
is 8-fold degenerate, because it has finite momentum, 
fe = (±27r/3,0),(0,±27r/3), and z-spin S z =±l/2; the 
single hole spectra of Figure 1 have been obtained by av- 
eraging over these 8 degenerate states. When energies 
are measured in units of J, it is immediately obvious 
that apart from an overall slight shift to lower energies 
in the 2-hole case (which may reflect the 'softening' of 
the spin excitations due to increased hole doping) the 
two spectra are remarkably similar: they show 'features' 
at comparable energies, with comparable spectral weight 
and with an anlogous dependence of the weight on J/t 
(analogous features can also be seen in the IPES spec- 
tra of the 16-site cluster). Up to the 'softening' of peak 
energies, the IPES spectrum thus scales with rih over its 
entire width, in strong contrast e.g. to noninteracting 
electrons, where the parts far from the Fermi energy re- 
main unaffected by a change of electron density. We thus 
see the same remarkable continuity with hole doping as 
established previously for the PES spectrum || and the 
dynamical charge correlation function J5| . 
We proceed to a detailed examination of the single hole 
case. More precisely, we study the addition of an f- 
spin electron to the single-hole GS with momentum 
(— 27r/3,0) and z-spin —1/2; the final states are eigen- 
states of the undoped antiferromagnet, i.e. spin ex- 
cited states, with z-spin 5 Z =0 and total spin St o t—0, 1. 
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We define the spin excitation operators Sx t j=Sj and 
S2,j=Sj ■ Sj+x+Sj ■ Sj+y (where e.g. j + x denotes the 
nearest neighbor of site j in x-direction) and their Fourier 
transforms S\^ q . Then, the following off-diagonal Green's 
functions describe the interference between IPES in the 
n + 1 hole GS and spin excitation of the n hole GS: 



B x (k,w)= J2 



(Qn\S{ tq \v> n }{Vn\c{jO n +l) 



^n\S{S X , q \Q n )^ 



5{u>-u>l). 



We have normalized the state Sx, q \O n ) to unity in or- 
der to suppress a possibly strong q dependence of its 
norm; in evaluating B\(k,uj) we moreover readjusted 
for each k the (arbitrary) relative phase between the 
ground states | Oq) and |0i) such that the frequency in- 
tegral of B\(k,Lu) is real and positive. The imaginary 
part then vanishes identically. Figure 2a shows A + (k,Lu) 
and 1 3? B x (k,u>)\. In the IPES spectrum A + (k,u>) there 
is precisely one peak with excitation energy w„=0 at the 
position of the 'hole momentum', (27r/3,0). The corre- 
sponding final state is the GS of the undoped system, so 
that this peak obviously describes the inverse process of 
hole creation in the half-filled GS. In a single particle pic- 
ture, it would be the only peak expected. In the actual 
IPEs spectrum, however, the bulk of weight is carried by 
peaks with finite excitation energy and comparison with 
B\(k,a) establishes a hierarchy of increasing spin exci- 
tation: the lowest group of the finite energy peaks cor- 
responds to a single spin wave left behind in the cluster, 
the next highest group to a 'bimagnon' and most prob- 
ably even higher peaks correspond to an even stronger 
spin excitation left behind by the annihilated hole. As 
regards the Fermi surface, it is clear that only the low- 
est state in this hierarchy, corresponding to 'complete 
annihilation' of the quasiparticle should be considered. 
The IPES spectrum thus consists of the 'hole pocket' 
at (27r/3, 0) H plus the 'shake off bands' at higher en- 
ergies. The latter are present throughout the Brillouin 
zone, however with a strongly reduced spectral weight 
near T; this is to be expected from the kinetic energy 
sum rule E kin = — J2k e kW£ PES , where denotes the 
noninteracting kinetic energy and W£ PES the integrated 
IPES weight for momentum k: large (small) IPES weight 
near the zone boundary (zone center) optimizes the ki- 
netic energy. 

We proceed to the two-hole case. The GS in this case 
is a spin singlet and has d x 2_ y i symmetry; introducing 
fci = (27r/3, 0) and k^ = (0, 2tv/3), the simplest two- 
particle state with these quantum numbers to be ex- 
pected within a rigid-band picture would be: 



-ka,l 



\vac) 



(1) 



For this state, the 'quasiparticle occupation' of the 4 
equivalent momenta ±fci, ±&2 is the same for each spin 



direction. Thus if we make the simplest assumption pos- 
sible, namely that the annihilation of one hole is not at 
all influenced by the presence of the second one, the two- 
hole IPES spectrum should simply be the average of the 
single hole IPES spectra for the different momenta and 
spin directions, as is indeed the case. Figure 2b shows 
the 'interference spectrum' \dlB\(k,uj)\ for n—1 (the lat- 
ter spectrum being averaged over the 4 degenerate single 
hole ground states with spin 1/2). Unlike the n=0 case, 
the IPES peaks now show overlap with both types of 
spin excitations - this is to be expected, because the ex- 
tra spin defects generated by the operators S\. q in the 
single hole GS may also be 'absorbed' into the dressing 
cloud of the hole. For example, addition of a spin ex- 
citation with momentum (27r/3, 27r/3) to the single hole 
GS with momentum (— 27r/3,0) with a certain probabil- 
ity may simply give the degenerate GS with momentum 
(0,27r/3). Apart from this slight complication, however, 
the structure of B\{k,uj) is quite similar to the single 
hole case: the single spin wave operator S\ has strongest 
overlap with the dominant peak at the bottom of the re- 
spective IPES spectrum, whereas the smaller IPES peaks 
at higher energies rather correspond to the 'bimagnon': 
obviously there is the same hierarchy of increasing spin 
excitation as in the single hole case (as is to be expected 
from the simple scaling of the IPES spectrum with rih, 
compare Figure 1. 

Taken together, the numerical data suggest the existence 
of different 'channels' in the IPES spectrum: in addi- 
tion to a 'conventional channel' at the lowest excitation 
energies, where a hole-like quasiparticle is 'annihilated 
completely', there exist higher energy final states where 
the annihilation process of the quasiparticle leaves be- 
hind a variable number of spin excitations, which formed 
the dressing cloud of the annihilated hole. The high 
degree of continuity of the IPES spectrum with doping 
moreover suggests that the spin-bag like quasiparticles 
for the half-filled case persist as well-defined entities also 
in the two-hole case (i.e. at nominal hole concentrations 
6~12%). Assuming that these different channels remain 
'disconnected' in the infinite system, i.e. that they do not 
merge to form a single 'band', the Fermi surface there- 
fore should take the form of hole pockets, consistent with 
increasing numerical evidence ||. The emerging picture 
of the low energy weight distribution and Fermiology in 
the full single particle spectral function then is summa- 
rized in Figure 3a (thereby the incoherent continua deep 
below the Fermi energy are omitted). To begin with, the 
Fermi surface is a hole pocket || generated by shifting 
the chemical potential into the more or less rigid next- 
nearest neighbor hopping band for hole motion in the 
undoped antiferromagnet Q). Thereby the PES spectral 
weight of the quasiparticle band in the outer parts of the 
Brillouin zone is small 0J(|; this 'shadow band' effect 
is readily understood B by the corollary of the kinetic 
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energy sum-rule, E kin =J2k £kW£ ES ', with W^ ES the in- 
tegrated PES weight. Low weight in the shadow band 
allows on one hand to maintain the next-nearest neigh- 
bor dispersion of width J, which seems to be optimal for 
the exchange mediated propagation of a single hole [ft), 
while at the same time staying close to the free-electron 
W£ ES (see Figure 3a) so as to optimize the 'ordinary' 
kinetic energy. We turn to the IPES spectrum, which 
consists of several components: the hole pockets right 
at Ep near (ir/2,ir/2) (or (n, 0) for two holes in a small 
clusters [||) and the various 'shake-off bands' at higher 
energy and predominantly in the outer part of the Bril- 
louin zone. When viewed with coarse energy and/or k 
resolution, the resulting distribution of spectral weight 
may be rather similar to a 'renormalized free electron 
band' jij; let us note that this is largely necessitated by 
elementary rules || . Key features which would allow for 
a distinction are a) the low intensity 'shadow bands' in 
the outer parts of the Brillouin zone, b) a 'disconnected' 
IPES spectrum consisting of low energy hole pockets near 
(-7r/2,7r/2) and higher energy magnon bands near (n 7 ir), 
and c) a Fermi surface with a volume proportional to the 
hole concentration 6. 

To address these issues, we proceed to a comparison with 
experiments on cuprate superconductors. PES exper- 
iments on the insulating antiferromagnet SriCuOiCl2 
have given evidence for a dramatic shadow band ef- 
fect in this undoped antiferromagnet. Assuming that the 
weight of the shadow bands even decreases in the doped 
case (as is suggested by cluster studies [§,[§) would read- 
ily explain the non-observation of the shadow bands in 
early ARPES experiments. A recent Fermi surface map- 
ping of the metallic compound -B^S^C 'aCu^Os+z with 
an extremely high density of fc-points by Aebi et al. [ jTcj ] 
has indeed shown indications for the shadow bands also 
in the doped case, a result which has been corroborated 
by conventional ARPES in the meantime . 
We turn to IPES, where experimental data are more 
sparse than in the PES case. IPES experiments which 
revealed several dispersive bands were reported by Bern- 
hoff et al. Jl^j . Figure 3b shows the positions of IPES 
peaks found by these authors in Bi2Sr2CaCu20%, as well 
as an assignment of bands compatible with our scenario. 
Grouping the higher lying peaks as indicated in the fig- 
ure, one can obtain almost quantitative agreement with 
the LDA band structure provided the two C11O2- 
plane derived bands are omitted from the latter. We are 
then left with only two 'disconnected' low energy fea- 
tures: the extreme low energy intensity near (7r/2,7r/2), 
which we associate with the hole pocket plus some 'iso- 
lated' spectral weight near (tt,tt). The latter feature is 
most clearly visible at the largest possible distance from 
fcf and seems to disappear when approaching kp. This 
is in rather clear contrast to the behaviour expected for 
a 'renormalized free-electron band', which on the con- 
trary should be most clearly visible near kp. On the 



other hand, this is precisely the behaviour expected for 
the magnon bands in the t—J model, whose large inten- 
sity near (tt, n) is favourable to lower the kinetic energy. 
We turn to the volume of the Fermi surface. Here the 
most direct probe the deHaas-vanAlphen effect, which is 
insensitive to the quasiparticlc weight Z and hence may 
detect the 'shadow' sheets of the Fermi surface. Mod- 
elling the holes as spin 1/2 particles distributed evenly 
over 4 symmetric pockets of equal size, the fraction of the 
Brillouin zone covered by a single pocket is 6/ 8. With a 
lattice constant of 3. 85 A (as would be appropriate for 
most cuprate superconductors) 5=0.18 gives a de Haas- 
van Alphen cross-section of 0.63/cT, which is roughly 
half-way in between two orbits of (0.53 ± 0.02)fcT and 
(0.78 ±0.02)fcT actually observed in YBa 2 Cu 3 6 . 97 @. 
This suggests to assume that the two 't—J bands' derived 
from the two CuC>2 planes/unit cell form bonding and an- 
tibonding combinations so that there is a disproportion- 
ation of holes between these two bands. With <5i=0.15, 
$2=0. 22 (i.e. an average 8 of 0.185) one obtains orbits 
of 0.52fcT and 0.76fcT, consistent with experiment fu) . 
For completeness we note that the experimental doping 
dependence of the low temperature Hall constant, which 
is another Z-independent quantity directly related to the 
carrier concentration, can also be reproduced well within 
the hole-pocket picture [ p^|JTc| 

. In summary, we have presented a cluster diagonaliza- 
tion study of the inverse photoemission spectrum for the 
low doping regime of the t—J model. The results first 
of all show a remarkable continuity with doping which, 
together with similar results for the photoemission spec- 
tra pj and dynamical correlation functions suggests 
that the single hole is the key problem for understand- 
ing the moderately doped region. The emerging picture 
of the single particle spectral function is that of a spec- 
tral weight distribution which in the neighborhood of the 
Fermi surface roughly resembles that of free electrons, 
combined with a hole-pocket Fermi surface with a vol- 
ume proportional to the number of doped holes. Various 
experimental data which suggest that this in cuprate su- 
perconductors are consistent with this scenario. 
Financial support of R. E. through the European Com- 
munity is most gratefully acknowledged. Computations 
were carried out at the Institute for Molecular Science, 
Okazaki. 

Figure 1: IPES spectrum A + (k 1 u})/rih for different k 
and J ft for n^=l (a) and n^=2 (b). ^-functions are re- 
placed by Lorentzians of width 0.1 J. 

Figure 2: (a) IPES spectrum A + (k,Lu) (full line), 
"interference spectra" 0.5 • Bi(k,uj) (dashed line) and 
B 2 {k,bj) (dotted line), for n=l. The ratio J/t=0.4, 5- 
functions are replaced by Lorentzians of width 0.01J. (b) 
Same as (a) for n=2. 
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Figure 3: (a) Schematic spectral function of the t— J 
model near the Fermi energy. The full line denotes the 
PES intensity from the quasiparticle band, the dashed 
line the IPES spectrum from the 'hole pocket', the 
short-dashed line the 'shake-off bands', (b) Compari- 
son of experimental and theoretical band structure of 
Bi^Sr^CaCu^O^,. The Fermi level is the zero of energy, 
the symbols denote prominent peaks in the IPES spec- 
trum (l2) (above E F ) and in the PES spectrum (be- 
low Ep). Dashed lines show the unoccupied part of the 
LDA bandstructure with the Cu02-plane derived 
bands removed, and rescaled and shifted according to 
E=0.77-(E LDA -E F )-Q.32eV. The full line indicates the 
H—J band'. A group of PES peaks originating from the 
superlattice structure of -B^SV^CaC^Os 11 1 is omit- 
ted. 
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